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LES  OVERVIEW 
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Dip.  di  Ing.  Aeronautica  e Spaziale 
Politecnico  di  Torino,  Italy 


Abstract.  In  the  paper  the  different  current  approaches  to  the  Large  Eddy 
Simulations  of  turbulent  flows  are  examined  and  critically  discussed.  A LES 
is  usually  defined  as  a numerical  computation  in  which  some  scales  are  cap- 
tured and  others  are  modelled,  but  different  procedures  are  used  for  such 
decomposition.  As  regards  the  modeling  of  the  unresolved  scales  they  are 
in  the  paper  tentatively  classified  as  exact  modeling  and  statistical  model- 
ing. A third  procedure,  the  dynamic  modeling  is  examined  as  a technique 
based  on  the  use  of  general  relations  that  connect  the  subgrid  scale  turbu- 
lent stresses  at  different  levels.  Their  implementation  to  standard  models, 
like  the  Smagorinsky  eddy  viscosity  one,  can  improve  the  results  and  can 
efficiently  extend  the  range  of  their  application.  The  connections  between 
LES  and  RANS,  the  simulations  based  on  the  Reynolds  Averaged  Navier- 
Stokes  equations,  are  finally  examined  and  some  critical  issues  related  to 
the  extraction  of  the  statistical  averages  from  LES  databases  are  discussed. 


1.  Introduction 

An  overview  of  the  many  approaches  that  can  be  classified  as  Large  Eddy 
Simulations,  (LES),  of  turbulent  flows,  cannot  be  done  without  considering 
the  limiting  procedures  of  the  Direct  Numerical  Simulation,  (DNS),  and 
the  Reynolds  Averaged  Navier-Stokes  Simulation,  (RANS).  If  compared 
with  the  complex  scenario  of  different  definitions,  filtering  operators,  trun- 
cations and  models  adopted  by  the  current  different  LES  techniques,  these 
two  asymptotic  approaches  to  the  simulation  of  turbulent  flows  should  ap- 
pear well  defined,  consolidated  basic  points,  but  that  is  far  from  true.  As 
regards  the  Reynolds  averages  the  only  thing  that  we  can  say,  following 
Frisch  (1995),  is  that  for  the  m,oment  the  partial  understanding  of  chaos  in 
deterministic  systems  gives  us  the  confidence  that  a probabilistic  descrip- 
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tion  of  turbulence  is  justified,  and  as  regards  the  DNS  it  is  important  to 
remark  with  Lee  (1995)  that  an  instantaneous  flowfield  of  direct  numerical 
simulation  must  be  viewed  not  as  the  true  solution  of  the  Navier  Stokes 
equations,  but  only  as  a snapshot  of  the  evolving  flowfield  at  some  fictitious 
time.  Uncomfortable  as  it  can  be,  anyway  we  can  only  conclude  with  Lee 
(1995)  that,  due  to  the  chaotic  nature  of  the  turbulent  field,  the  only  viable 
outcome  of  both  approaches,  DNS  and  RANS,  is  prediction  of  the  aver- 
aged flow  quantities,  and  that  is  only  the  beginning  of  a lot  of  problems.  A 
statistical  description  needs  the  definition  of  a statistical  ensemble  of  flow 
realizations,  and  following  with  the  citations  let  us  notice  now  with  Aubry 
(1991)  that  this  definition  in  the  applications  has  been  quite  flexible.  As  a 
matter  of  fact  the  statistical  ensemble  usually  consists  of  a temporal  domain 
under  the  ergodicity  assumption,  or  of  a symmetry  group  under  which  the 
equations  for  that  particular  flow  are  invariant,  or  of  an  ensemble  of  dis- 
crete times  conditionally  determined.  The  basic  ensemble,  a set  of  initial 
conditions,  is  usually  never  considered,  for  obvious  practical  considerations. 

Let  us  now  examine  the  LES.  Prom  the  basic  point  of  view  a Large 
Eddy  Simulation  should  produce  a filtered  representation  (a)  / 

(a)f  = F(a)  (1) 

of  the  original  quantity  a,  where  T is  a general  filtering  operator  provided 
with  some  particular  smoothing  properties.  One  particular  aim  of  this  ap- 
proach is  to  derive,  given  the  evolutionary  equation  of  a,  the  evolutionary 
equation  for  the  filtered  quantity  (a)/.  That  is  not  so  easy,  and  different 
approaches  can  be  imagined.  Let  us  define  as  usual,  Leonard  (1974),  the 
filtered  quantity  (a) $ as  given  by 

(a)f  = J a(t')F(t  — t')dt'  (2) 

where  F{z)  is  a frequency  function,  see  Hirschman  & Widder  (1955).  If 

da  t / \ 

dt = A{o)  (3) 

is  the  evolutionary  equation  for  a,  the  evolutionary  equation  for  (a)  / 

will  usually  contain,  if  non  linear,  unresolved  moments.  Three  basic  ap- 
proaches can  be  conceived,  one  strictly  deterministic  and  the  others  sta- 
tistical in  order  to  close  this  filtered  equation.  The  first  one,  that  we  will 
call  exact  modeling,  is  based  on  the  analytic  properties  of  the  convolution 
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kernel,  the  frequency  function  F(z ),  and  it  is  clear  that  for  very  high  res- 
olutions a LES  is  similar  to  a DNS,  so  that  statistical  ingredients  seem 
unadapted.  We  remark  that  the  statistical  approaches  are  based  on  a prob- 
abilistic interpretation  of  the  definition  (2),  and  that  seems  a little  difficult 
to  understand.  It  is  however  easy  to  see  that  we  can  read  this  physical  de- 
terministic average  in  terms  of  a filtered  density  function  py(a,  t).  We  can 
write 

(■ a)/  — J pf(a,t)ada  (5) 


where 


pj(a,t)  = J 5[a  — a(t)]F(t  — t')dt' 


and  two  statistical  procedures  can  be  conceived,  the  first  one  based  on 
the  derivation  of  the  evolutionary  equations  for  the  different  moments  that 
appear  in  the  evolutionary  equation  of  (a)  f and  the  second  one  based  on 
the  evolutionary  equation  for  the  filtered  density  function  pj(a,t).  That  is 
similar  to  RANS.  Also  in  this  case  the  usual  definition  of  the  mean,  easy 
to  compare  with  the  experiments,  is  based  on  an  infinite  time  average,  but 
from  the  fundamental  point  of  view  it  is  more  correct  to  think  in  terms  of 
an  ensemble  operator  £ based  on  a set  of  initial  conditions 


(o)e  = J Pe(a)a(a,t) da 


where  Pe,(a)  is  a probability  density  function  on  the  initial  states  at  the 
time  t = 0 

a = a(0)  (8) 

If 

da  a / \ 

It  = A(a)  (9) 

is  the  evolutionary  equation  for  a,  also  in  this  case  the  basic  problem  is  to 
derive  the  evolutionary  equation  for  ( a.)e 


and  we  remark  that  if  we  write 


(a)e  = J pe{a,  t)ada 


where  now  pe(a , t ) is  a probability  density  function  on  a,  given  by 


pe(a,t)  = J S[a  - a(a,t.)]Pe(a)dc 
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another  challenging  problem  is  to  derive  the  evolutionary  equation  for 

Pe(a,t) 

dpe 

dt 

In  the  following  we  will  examine  these  two  approaches,  the  exact  model- 
ing and  the  statistical  modeling , in  some  detail  and  we  will  finally  make 
some  comments  on  a third  modeling  technique  that  is  of  some  usefulness 
in  improving  dynamically  the  previous  procedures.  We  will  call  it  dynamic 
modeling  and  we  remark  that  it  is  based  on  general  operational  relations 
that  any  modeling  procedure,  deterministic  or  not,  must  verify,  like  the 
tensorial  or  the  Galilean  invariance. 


2.  Exact  modeling 

We  tentatively  classify  as  exact  modeling  the  procedures  based  on  the  an- 
alytical properties  of  an  explicit  convolutional  filtering  operator  T 

/OO 

a(t')F(t  - t')dt'  (14) 

-OO 


where,  as  remarked  in  the  introduction,  the  kernel  F(z ) is  a frequency 
function.  It  is  interesting  to  notice,  Hirschman  & Widder  (1955),  that  if 
L(s)  is  the  inverse  of  the  bilateral  Laplace  transform  of  F(z) 


(15) 

a useful  inversion  formula  is  the  following 

L(T>)(a)  f = a 

(16) 

where  V stands  for  the  derivative  operator 

III 

9r|a. 

(17) 

From  the  operational  point  of  view  the  differential  operator  L(V)  is  the 
inverse  of  the  integral  operator  T 


L(V)  = F~l  (18) 

and  an  interesting  class  of  filters  is  the  following 

T~x  = (1  - ai©)((l  - a2V)  • • • (1  - anV ) (19) 

It  can  be  shown  that  these  filters  are  variation  diminishing  in  the  sense  that 
the  number  of  changes  of  sign  of  the  filtered  function  (a)  j never  exceeds  the 
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number  of  changes  of  sign  of  the  original  function  a,  and  other  interesting 
properties  are  related  to  the  mean  my  and  the  variance  Vj  of  the  kernel  of 
the  inverse  convolutional  filter,  the  frequency  function  F(z),  given  by 

(20) 

We  remark  that  the  smoothing  properties  of  the  variation  diminishing 
transform  are  particularly  interesting  for  LES,  and  till  now  unexplored. 
Historically  they  were  first  introduced  by  Schoenberg,  and  are  connected 
to  the  zero  crossing  rate  of  a turbulent  fluctuation  that,  as  remarked  in  the 
past  by  Liepmann,  is  directly  related  to  the  dissipation. 

Many  other  different  filters  have  been  proposed  for  LES,  and  a complete 
review  of  all  them  should  require  many  pages.  We  refer  for  that  to  the  recent 
book  of  Sagaut  (2001)  where  definitions  and  properties  of  classical  filters 
for  LES  are  discussed  both  for  the  homogeneous  and  the  inhomogeneuos 
case.  This  last  issue  is  very  important  for  the  applications,  and  we  remark 
that  a general  class  of  differential  filters,  Germano  (1986),  is  particularly 
interesting  in  order  to  explore  what  happens  when  nonhomogeneous  filter- 
ing operators  do  not  commute  with  differentiation,  Germano  (2000).  As  an 
example  let  us  consider  a linear  differential  form  of  the  second  order 


mj 


n 1 

E- 

1 Oik 


«/  = Ei 


i ak 


V — X + dVt  + A*,!?*;  + A kiDk^i  (21) 


where 


Vt  = 


dt 


Vk  = 7T 

dxk 


(22) 


and  let  us  associate  to  this  differential  operator  a filtering  operator  fF  de- 
fined as 

X = VT  (23) 


where  X is  the  identity.  It  is  interesting  to  remark  that  the  differential 
operator  V is  by  definition  the  inverse  of  the  filtering  operator  T 


V = T~l 


(24) 


and  formally  T can  be  expressed  by  an  integral  convolution  with  a ker- 
nel given  by  the  Green’s  function  associated  to  V.  We  notice  that  if  the 
parameters  d,Ak,Aki  are  not  constant  the  operator  T does  not  commute 
with  the  derivatives  T>t  and  Vk,  and  what  is  interesting  with  these  filters 
is  that  provided  with  the  exact  inverse  of  T we  can  easily  write  the  exact 
expressions  for  the  commutative  terms 


TVt  - VtT 


xv  k — VkX 


(25) 
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In  fact  if  Q is  a generic  operator  we  obtain  by  applying  the  relation  (23) 

QVT  = VTQ  (26) 

and  we  can  write 

(QV  - VQ)T  = V{TG  - GT)  (27) 

If  we  now  compare  this  relation  with  the  relation  (23)  we  obtain 

TG-QT  = T{QV  - VG)T  (28) 

that  provides  the  exact  form  for  the  commutative  errors  in  terms  of  the 
filtered  quantities. 

3.  Statistical  modeling 

As  remarked  before  statistical  modeling  reads  the  filtering  length  as  an 
interval  of  indeterminacy  that  separates  what  can  be  calculated  analytically 
from  what  has  to  be  guessed  statistically.  This  point  of  view  can  be  similarly 
applied  to  a projection.  If  we  express  the  generic  turbulent  velocity  field  Ui 
in  terms  of  a generalized  Fourier  expansion 

00 

ui  = Yl  U^k  (29) 

1 

where  tpk  is  a particular  set  of  basis  functions  and  are  random  Fourier 
coefficients,  following  Yoshizawa  (1982)  we  can  define  a partial  statistical 
operator  £/  that  applied  to  it*  gives  the  J7- level  statistical  representation 

/ 00 

£fiui)  = y ^ uikfk  + (uik)elPk  (30) 

1 /+1 

The  partial  statistical  theory  of  turbulence  is  at  the  beginning,  and  in 
the  opinion  of  the  author  the  Large  Eddy  Simulation  has  stimulated  the  re- 
search in  the  field.  We  recall  that  the  probability  density  function  approach 
has  been  pioneered  in  turbulence  by  Lundgren  (1967).  His  method  derives 
directly  the  evolutionary  equation  for  the  ensemble  pdf  from  the  differen- 
tial equations  which  define  the  conservation  laws,  and  the  principal  fields 
of  application  are  the  reacting  turbulent  flows,  but  recently  this  approach 
has  been  extended  to  LES,  see  Pope  (1990)  and  Gao  &;  O’Brien  (1993). 
This  concept  of  a pdf  within  the  subgrid , Madnia  &:  Givi  (1993),  is  very 
promising  for  LES  and  some  first  simulations  based  on  a velocity  filtered 
density  function  computed  by  a Lagrangian  Monte  Carlo  procedure  have 
been  recently  performed  by  Gicquel  et  al.  (2001). 
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As  remarked  in  the  introduction  the  exact  and  the  partial  statistical 
modeling  based  on  filtered  pdf  could  seem  incompatible.  If  the  filtering 
operation  is  defined  as  a convolution  in  the  physical  space 


(a)f  = J a(x')F(x  — x')dx 


(31) 


we  have 

Tf(a,b)  = Tf(—a,  —b)  (32) 

where  by  definition  the  generalized  central  moment  associated  to  the  quan- 
tities a and  b is  given  by 

Tf{a,b)  = {ab)f  - (a)f{b)f  (33) 

and  this  condition  of  reversibility  should  be  respected  by  modeling.  It  is 
however  important  to  recall  that  the  Smagorinsky  model,  based  on  sta- 
tistical considerations  and  applied  to  the  anisotropic  part  of  the  subgrid 
turbulent  stress, 


Tf(ui,Uj)  ~ Mf(ui,Uj ) = -2uf(sij)f 


(34) 


is  such  that 

M/(ui,Uj ) = - Mf(-~Ui,-Uj ) (35) 

and  this  condition  is  violated.  This  situation  is  reminiscent  of  the  reversibil- 
ity paradox  related  to  the  connections  between  the  statistical  kinetic  theory 
and  the  microscopic  laws  of  mechanics  that  are  invariant  under  time  rever- 
sal. We  remark,  Germano  (2001),  that  if  we  rewrite  the  integral  (31)  in  its 
Lebesgue  version 


where 


(a) f = J a(x')F(x  — x')dx'  = J ap(x , a)da 
p(x,a)  — J F(x  — x')6[a  - a(x')]dx' 


(36) 


(37) 


the  Smagorinsky  model  could  be  recovered  by  an  assumed  probability  dis- 
tribution p(x,  Ci)  given  by 


p{x,Ci)  = (27rt/2)-3/2exp[-C2/2I/2 


1-22. 

u 2 


C2  5 \ CidU2 + Ci{sij)fCj 


2U2  2)  U2  dx 

where  — ut  — (uz) /,  G'2  = CtCi  and 

3U2  = Tf(ui,Ui ) 


U2 


(38) 


(39) 
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It  seems  in  conclusion  that  the  eddy  viscosity  model  could  only  be  justified 
on  the  basis  of  probabilistic  arguments,  and  that  the  irreversible  LES  could 
only  find  its  formulation  in  the  large  eddy  probabilistic  density  function. 


4.  Dynamic  modeling 


The  dynamic  modeling  procedures  are  based  on  identities  that  relate  the 
subgrid  scale  stresses  at  different  resolution  levels.  The  simplest  one  is  the 
following 

Tfg{Ui,Uj)  = (TfiUUUj^g  + Tg((Ui)f,  (Uj)f)  (40) 

where  the  various  terms  are  given  by  definition  by  the  expressions 


T f (ui , Uj  ) 
Tfg{Ui,Uj) 
(Tf(Ui,Uj))g 
Tg{{Ui)f,  ( Uj)f  ) 


!F(uiUj)  — ^(u^Tluj) 

QT(uiUj)  - GT{ui)QT(uj) 

GT(uiUj)  - Gi^iu^Tiuj)) 
g(T{ui)T{Uj))  - gF(Ui)gF(Uj) 

(41) 


and  where  T and  G are  respectively  the  LES  filter  and  the  test  filter.  Also 
in  this  case  to  discuss  the  many  applications  and  variants  adopted  since  the 
first  dynamic  model,  Germano  et  al.  (1991),  should  require  a lot  of  space, 
and  we  refer  to  a recent  paper  of  Piomelli  (1999)  for  a critical  analysis. 
It  is  perhaps  more  interesting  for  this  brief  review  to  remark  that  this 
procedure  is  based  on  purely  formal  properties  and  can  only  improve  an 
existing  model,  like  the  Smagorinsky  eddy  viscosity  model,  or  the  structure 
function  model.  That  is  good  from  one  side,  because  as  such  its  range  of 
application  is  very  large,  and  we  recall  a recent  interesting  extension,  Im 
et  al.  (1997),  of  the  dynamic  procedure  applied  to  the  so  called  G-equation 
that  predicts  the  evolution  of  the  propagating  front  of  laminar  flamelets 
corrugated  by  turbulent  eddies.  From  the  other  side  however  that  is  a limit, 
because  in  some  sense  there  is  nothing  new  that  is  added  physically  by  this 
procedure,  and  the  improvement  is  due  to  a better  consistency  with  the 
formal  properties  of  the  real  subgrid  stresses.  This  point  seems  to  the  author 
very  important  and  has  been  the  matter  of  an  interesting  criticism  by  Pruett 
(1997).  With  reference  to  the  relation  (40)  he  states  that  this  identity  is 
mathematically  tautological,  and  unnecessary  as  a basis  for  dynamic  SGS 
models , and  in  our  opinion  this  criticism  is  applicable  or  not  according  to 
the  situation.  Let  us  examine  the  identity  (40)  from  the  point  of  view  of 
the  exact  modeling.  In  this  case  we  assume  that  both  the  test  filter  G and 
the  LES  filter  T are  explicitly  known,  and  as  in  Brun  & Friedrich  (2000) 
we  can  expand  them  in  one  dimension  as  follows 

J~  — X + QfV  + bjT)^  + CjTft  • • • 
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Q 'L  + dg'D  + bgV2  + CgT)^  • • • (42) 

where  V stands  for  the  derivative.  If  the  problem  is  to  find  an  expression 
for  the  turbulent  stress  at  the  TQ 


Tfg(ui,uj)  = GFiuiUj)  - GF(ui)GT(uj ) (43) 

we  can  in  this  case  derive  an  explicit  expression  for  TQ 

GT  = l + afgV  + bfgV2  + cfgV 3---  (44) 


where 


af  + ag  = 0,fg 

bf  + dfdg  + bg  = bfg 

Cf  + bjdg  + dfbg  +Cg  - Cfg 

(45) 


and  it  is  clear  that  in  this  case  Pruett  is  right.  In  fact  we  can  easily  derive 
directly  that 


TfgiuuUj)  = {2b fg  - a2fg)V(ui)T>(uj)  + ■ ■ • (46) 

and  the  identity  (40)  is  unessential.  On  the  contrary  if  we  are  provided 
only  with  an  explicit  LES  model , and  we  do  not  exactly  know  the  explicit 
originating  LES  filter,  we  can  only  write  that 


dTf{ui,Uj)  , L d2Tf{ui,uj ) 

(Tf(lli,Uj))g  — Tf(Ut,Uj)  + dg  +bg 

Tg((ui)f,{uj)f)  = (2bg  - a2)V{{ui) f)V{{uj) f)  + ••• 


+ 


(47) 


and  the  only  way  to  derive  the  turbulent  stress  at  the  TG  level  remains  to 
apply  the  identity  (40). 


5.  Connections  between  LES  and  RANS 

The  comparisons  among  DNS,  LES  and  RANS  databases  pose  a lot  of 
problems  that  in  some  cases  can  be  neglected,  but  that  conceptually  remain. 
The  first  concerns  the  recovery  of  the  statistical  data  from  a LES  database. 
We  remark  that  what  we  produce  by  LES  are  filtered  quantities,  for  example 
the  filtered  values  of  the  velocity  field,  (tq) /.  When  we  apply  an  explicit 
model  we  have  no  idea  of  its  explicit  generating  filter,  if  ever  exist.  If  we 
define  the  LES  fluctuations  u[  to  the  statistical  mean  as 

ui  = (ui)e  ~ iui)f 


(48) 
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we  notice  that  only  in  the  case  (u'i)e  = 0 the  simple  usual  relations 

Te(Ui,Uj)  = {u'iu'j)e  + {Tf(Ui,Uj))e 
Te{Ui,Uj,Uk)  = (u'iu'ju'k)e  + {Tf{Ui,Uj,Uk))e  + (UiTf(uj,Uk))e  + 

+ (UjTf{uk,Ui))e  + (u'kTf(Ui,Uj))e  (49) 

are  valid,  and  that  is  not  so  good.  As  remarked  in  Hussaini  et  al.  (1989) 
if  higher  order  turbulence  statistics  are  needed  beyond  the  mean  velocity, 
the  problem  of  defiltering  arises,  and  this  problem  becomes  most  critical 
when  more  than  10  — 20%  of  the  turbulent  kinetic  energy  is  in  the  subgrid 
scale  motions.  The  literature  on  the  estimation  of  the  second  and  the  third 
order  statistical  moments  based  on  LES  databases  is  relatively  poor,  and 
in  some  cases  we  are  obliged  to  think  that  probably  the  subgrid  models 
are  considered  very  important  dynamically,  in  order  to  compute  the  mean 
values,  but  not  so  reliable  in  order  to  be  implemented  in  the  relations  (49). 
We  remark,  Germano  (2001),  that  only  in  the  case  that  the  following  chain 
of  relations  is  satisfied 


«}e  = 0 

(Tf(Ui,Uj))e  = 0 

(uir f (uj > uk)) e = 0 

(r/(Ui,UJ-,Ufc))e  = 0 


(50) 


the  filtering  procedure  produces  statistical  moments  simply  given  by 

Te(Ui,  Uj)  — (UjUj) e 
Te(ui,Uj,Uk)  = (u^UjUk)e 
Te(Ui,Uj,Uk,Ul)  = (u-u'UfcU|)e 

...  = •••  (51) 

and  the  study  of  explicit  filters  or  explicit  models  that  implement  these 
statistical  constraints  are  currently  under  way. 

6.  Conclusions 

The  various  modeling  approaches  and  numerical  procedures  that  since  the 
first  applications  in  meteorology  have  been  applied  to  simulate  the  large 
scales  of  turbulence,  have  stimulated  during  the  time  a lot  of  different  in- 
terests and  have  raised  a lot  of  different  questions.  A general  overview  of  all 
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that  is  nearly  impossible,  and  inevitably  is  a little  biased  by  the  particular 
attitudes  or  tendencies  of  the  reviewer.  If  we  look  at  the  agenda  of  this 
meeting  we  can  see  that  the  applications  of  LES  are  very  promising  and 
wide-spread  on  a lot  of  different  fields  of  application,  from  the  aeronautical 
engineering  to  the  estimate  of  the  pollution  dispersion,  in  reactive  flows  as  in 
heat  transfer.  Industrial  needs  and  interests  motivate  obviously  the  future 
of  LES,  that  will  probably  still  remain  for  a long  time  the  only  viable  way 
for  accurate  calculations  of  turbulent  flows.  Hybrid  RANS-LES  strategies 
like  the  Detached  Eddy  Simulation  or  time  dependent  RANS  are  presently 
developed  in  order  to  simulate  complex  turbulent  flows  and  to  combine 
LES  with  standard  and  well  developed  computational  and  modeling  proce- 
dures. In  this  particular  case,  as  in  many  others,  practical  and  theoretical 
issues  are  strongly  intermixed.  From  one  side  the  motivation  of  these  hy- 
brid procedures  is  mainly  practical  : LES  remains  after  all  very  expensive, 
and  simplest  computational  strategies  remain  very  attractive.  From  another 
side  however  the  final  goal  is  very  ambitious  and  aims  to  construct  a unified 
turbulence  model  useful  for  all  purposes,  from  wall  modeling  to  the  wake. 
This  point,  as  many  others,  is  a real  challenge  and  imposes  a big  effort  in 
a lot  of  different  fields,  from  the  exploration  of  higher  order  computational 
schemes  to  the  basic  research  on  turbulence,  from  the  study  of  unstruc- 
tured grids  to  the  speculations  about  statistical  and  deterministic  issues  in 
chaotic  systems.  The  fundamentals  of  turbulence  remain  one  of  the  most 
fascinating  unresolved  problem  and  the  hopes  placed  in  a more  rational 
design  of  aircraft,  automobiles,  turbines  and  many  others  industrial  pro- 
cesses based  on  an  improved  computational  ability  of  fluid  flows  are  more 
and  more  increasing.  In  order  to  cope  with  all  these  problems  we  need  the 
effort  of  everybody  and  a strong  cooperation  among  different  partners. 
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